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COMPRESSION CREEP OF FILAMENTARY COMPOSITES 

ABSTRACT 

Axial and transverse strain fields induced in 
composite laminates subjected to compressive creep loading 
were compared for several types of laminate layups. Unidirectional 
graphite/epoxy as well as multi-directional graphite/epoxy and 
graphite/PEEK layups were studied. Specimens with and without 
holes were tested. The specimens were subjected to compressive 
creep loading for a 10 hour period. All tests were conducted at room 
temperature and ambient humidities. 

In-plane displacements were measured using moire interferometry. 
A computer based data reduction scheme was developed which 
reduces the whole-field displacement fields obtained using moire to 
whole-field strain contour maps. This approach proved to be an ideal 
method of displaying the data obtained using moire interferometry. 
Only slight viscoelastic response was observed in matrix-dominated 
laminates, except for one test in which catastrophic specimen 
failure occurred after a 16-hr period. In this latter case the 
specimen response was a complex combination of both viscoelastic 
and fracture mechanisms. No viscoelastic effects were observed for 
fiber-dominated laminates over the 10 hour creep time used in this 
study. The experimental results for specimens with holes were 
compared with results obtained using a finite-element analysis. The 



comparison between experiment and theory was generally good. 
Overall strain distributions were very well predicted. The finite 
element analysis typically predicted slightly higher strains values 
at the edge of the hole, and slightly lower strains values at 
positions removed from the hole, than were observed experimentally. 
It is hypothesized that these discrepancies are due to nonlinear 
material behavior at the hole edge, which were not accounted for 
during the finite-element analysis. 


TABLE OF CONTENTS 


Page 

ABSTRACT ii 

LIST OF FIGURES v 

LIST OF TABLES vii 

ACKNOWLEDGEMENTS viii 

INTRODUCTION 1 

OBJECTIVES 4 

TESTING EQUIPMENT 6 

Compression Fixture 6 

Laser and Optics 12 

Camera and Film 13 

Specimen Preparation 1 3 

Test Matrix 14 

Selection of Creep Stress Levels 15 

CALCULATION OF STRAIN FROM MOIRE FRINGES 1 9 

Data Reduction by Mechanical Differentiation 19 

Data Reduction by the Displacement-Field Approach 21 

Verification of the Data Reduction System 32 

RESULTS 36 

SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 77 

Summary 77 

Conclusions 78 

REFERENCES 85 



APPENDIX A 

COMPUTER PROGRAMS 
Program MOIRE 
Program STRREG 


LIST OF FIGURES 


Number Page 

1. Compression Fixture in the SATEC Creep Frame. 

(Schematic) 8 

2. Mirror Arrangement for the Moire Fixture. (Schematic) 9 

3. Moire and Compressive Fixture in the SATEC Creep Frame 1 0 

4. Moire Optical Set-up 1 1 

5. Displacement-Position Plots For 3u/3x and 3u/3y 22 

6. Displacement-Position Plots For 3v/3x and 3v/3y 23 

7. V-Displacement Field, Aluminum Specimen. (Original) 25 

8. V-Displacement Field, Aluminum Specimen. (Cleaned) 26 

9. Fringe Center Locations 29 

10. Axial Strain Field Map, Aluminum Specimen 34 

11. Axial Strain Field Maps, Infinite Plate, Finite-Element 3 5 

12. V-Displacement Field, [0]48 t - 1 min 44 

13. U-Displacement Field, [0]48 t - 1 min 45 

14. Axial Strain Field Map, [0]48 t= 1 min 46 

15. Axial Strain Field Map, [0]48 t- 10 hrs 46 

16. Transverse Strain Field Map, [0]48 t* 1 min 47 

17. Transverse Strain Field Map, [0]48 t - 10 hrs 47 

18. Axial F. E. Strain Field Map, [0]48 48 

19. Transverse F. E. Strain Field Map, [0]48 48 

20. Axial Strain Field Map, [90]48 t = 1 min 52 

21 . Axial Strain Field Map, [90]48 t - 10 hrs 52 

22. Transverse Strain Field Map, [90]48 t - 1 min 53 

23. Transverse Strain Field Map, [90]48 t ■ 10 hrs 53 

24. Axial F. E. Strain Field Map, [90]48 54 

25. Transverse F. E. Strain Field Map, [90]48 54 

26. V-Displacement Field, [0/45/0/-45]6s t = 1 min 58 

27. U-Displacement Field, [0/45/0/-45]6s t * 1 min 59 

28. Axial Strain Field Map, [0/45/0/-45]6s t = 1 min ...60 

29. Axial Strain Field Map, [0/45/0/-45]6s t - 10 hrs 60 


30. Transverse Strain Field Map, [0/45/0/-45]es t = 1 min 61 

31. Transverse Strain Field Map, [0/45/0/-45]6 S t- 10 hrs 61 

32. Axial F. E. Strain Field Map, [0/45/0/-45] 6s 62 

33. Transverse F. E. Strain Field Map, [0/45/0/-45]6 S 62 

34. Delamination of Type D [90/-45/90/45]6s Gr/Ep 

V-Displacement 63 

35. Delamination of Type D [90/-45/90/45]6s Gr/Ep 

U-Displacement 64 

36. Front View of Failed [90/-45/90/45]6s Gr/Ep Specimen 65 

37. Side View of Failed [90/-45/90/45]6s Gr/Ep Specimen 65 

38. Axial Strain Field Map, [90/-45/90/45]6s t - 1 min 66 

39. Axial Strain Field Map, [90/-45/90/45] 6s t = 10 hrs 66 

40. Transverse Strain Field Map, [90/-45/90/45]6s t = 1 min 67 

41. Transverse Strain Field Map, [90/-45/90/45]6s t = 10 

hrs 67 

42. Axial F. E. Strain Field Map, [90/-45/90/45] 6s 68 

43. Transverse F. E. Strain Field Map, [90/-45/90/45]6s 68 

44. Axial Strain Field Map, 1/16" Dia Hole, Gr/PEEK, t = 1 

min 71 

45. Axial Strain Field Map, 1/8" Dia Hole, Gr/PEEK, t = 1 min 71 

46. Axial Strain Field Map, 1/4" Dia Hole, Gr/PEEK, t = 1 min 72 

47. Axial Strain Field Map, 3/8" Dia Hole, Gr/PEEK, t = 1 min 72 

48. Transverse Strain Field Map, 1/16" Dia Hole, Gr/PEEK, t = 1 

min 73 

49. Transverse Strain Field Map, 1/8" Dia Hole, Gr/PEEK, t = 1 

min 73 

50. Transverse Strain Field Map, 1/4" Dia Hole, Gr/PEEK, t = 1 

min 74 

51. Transverse Strain Field Map, 3/8" Dia Hole, Gr/PEEK, t = 1 

min 74 

52. Axial F. E. Strain Field Map, 1/4" Dia Hole, Gr/PEEK 75 

53. Transverse F. E. Strain Field Map, 1/4" Dia Hole, Gr/PEEK 75 


LIST OF TABLES 


Number Page 

1. Ultimate Strength 17 

2. Test Matrix 18 

3. Measured and Predicted Moduli for Gr/Ep Laminates 41 


ACKNOWLEDGEMENTS 


The authors wish to thank the many individuals who have helped 
make this work possible. Financial support was provided through 
NASA-Langley Grant No. NAG-1 -621. Mr. John D. Whitcomb has acted 
as Project Monitor. Sincere thanks are extended to Mr. Whitcomb for 
his encouragement and many helpful suggestions. The technical 
advice and insight of Prof. Ken Y. Lin is gratefully acknowledged. 
Thanks are extended to Profs. Minoru Taya and Colin H. Daly who 
edited the draft version of this report. Thanks is given to 
Mr. Raymond J. Klein for his technical advice regarding the use of 
the moire interferometry and loading system. Mr. Russ McGuire of 
Boeing Commercial Airplane Company provided the Gr/Ep panels used 
in the study. Profs. Dan Post and Robert Czarnek of Virginia Tech, 
Prof. Jon Epstein of Georgia Tech, and Dr. Vance Deason of Idaho 
National Laboratories all made valuable contributions to this study. 
And finally sincerest thanks goes to Mr. Mahyar S. Dadkhah and 
Visiting Prof. F.X. Wang of UW for their technical advice and 
assistance on all parts of this project. 



INTRODUCTION 


The advent of fiber-reinforced composites materials has been 
termed the "biggest technical revolution since the jet engine."! 1 1 
Advanced polymeric composite materials exhibit high 

strength-to-weight and stiffness-to-weight ratios. These qualities 
make composite materials extremely useful in applications where 
high performance and weight savings are important. Aircraft and 
spacecraft are typical weight-sensitive structures in which 
composite materials are cost effective. Composite materials 
usually exhibit the best qualities of their constituent components 
and often some desirable qualities that neither component possesses 
separately. Modern composites offer a dramatic improvement in 

mechanical/thermal capabilities, when compared to conventional 
monolithic materials, including an increase in; strength, stiffness, 
corrosion resistance, wear resistance, attractiveness, weight, 
fatigue life, thermal conductivity and acoustical insulation.! 1 1 

Composites made from strong continuous fibers embedded in a 
relatively weak matrix material exhibit anisotropic behavior when 
subjected to loading. Hence the mechanical and thermal behavior of 
composites is often much different than the behavior exhibited by 
familiar isotropic materials such as steel or aluminum. A 
considerable study of composites has focused on the response to 
tensile loading. Significantly less study has been devoted to the 
behavior of composites in compression. The research described 
herein was an investigation of the response of composite laminates 
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to uniaxial compressive loading, and should contribute towards a 
greater understanding of the compressive behavior of composites. 

Polymeric composites exhibit viscoelastic behavior which is 
more pronounced than that exhibited by traditional engineering 
materials. These effects, in general, become more evident at 
elevated temperatures, although significant viscoelastic behavior 
also occurs at room temperatures over long time periods. Exposure 
to various plasticizers such as moisture, jet fuel, or gasoline can 
also result in an increase in viscoelastic behavior. The expected 
service life of composite structures can be on the order of tens of 
years, and therefore viscoelastic response may result in delayed 
failure or excessive deflections after some years in service. 
Consequently a clear understanding of the nature and magnitude of 
this viscoelastic behavior is important with the increasing use of 
composites as substitutes for traditional materials. 

A study of composites in compression was initiated at the 
University of Washington in 1986. The study was supported through 
a two-year grant provided by the Fatigue and Fracture Branch of the 
NASA-Langley Research Center. The program can be roughly divided 
into two complementary efforts, one involving an analytic (i.e., 
finite-element) analysis and the second involving an experimental 
analysis. The present report describes the experimental efforts 
which have been performed during the second year of the 
UW/NASA-Langley study. Comparisons between experimental 
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results the finite-element analyses are presented where 
appropriate. 

The experimental portion of the study has involved the 
application of compressive creep loads to composite laminates. All 
tests have been performed at room temperatures and ambient 
humidity conditions. Two distinct material systems have been 
investigated: Hercules graphite-epoxy (G/E) IM7-8551 and Imperial 
Chemical Industries (ICI) graphite-PEEK (G/P) APC-2. Several 
different laminate layups and specimen geometries have been used. 

One unique aspect of the experimental analyses has been the 
use of moire interferometry to measure in-plane surface 
displacements. Moir6 interferometry is a state-of-the-art optical 
method which provides whole-field displacement contours. The 
technique was developed during the late 1970s' and early 
1 980s’f2]C3], Moire interferometry has since been used extensively 
to study the behavior of compositesl 4 H 5 H 6 ] as well as other 
structural materials( 7 H8][9][iO][i i]. However, to the authors' 
knowledge this is the first time moire interferometry has been used 
under creep conditions. 



OBJECTIVES 


The experimental efforts performed during the first year of 
the study have been previously described in several publications by 
Klein and Tuttle! 1 2 H' | 3][i 4] Klein's study was devoted primarily to 
the design and development of the experimental apparatus required. 
Specifically, Klein developed a compressive loading fixture (required 
to apply a uniaxial compressive creep load to the composite 
specimen) and a laser-based optical system (used to measure the 
resulting surface displacements using moire interferometry). The 
experimental system was successfully used by Klein to obtain moire 
patterns for a relatively limited number of tests. In-plane strains 
were then calculated using an essentially graphical procedure. This 
technique was tedious and time-consuming. Furthermore, the data 
reduction — 3ohome — e - nty - allowed for calculation — of — strains — along 
discrete lines or "slices" through the moire image! 12 !. Hence, 

whole-field strain measurements were not obtained. Klein 
recommended that an alternate data reduction scheme be developed 
to automate the data reduction process, hopefully resulting in a 
more accurate calculation of strain levels as well as whole-field 
strain contours. 

The present report describes the experimental efforts 
performed during the second year of the UW/NASA-Langley study. An 
expanded series of tests have been performed, and a more 
sophisticated data reduction scheme has been developed. The 
specific goals of the present study were as follows: 
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• To perform creep tests using unidirectional [0]4 8 
and [90]48 laminates. These data would allow 
calculation of the engineering moduli E-n, E22, n-12, 
and would also indicate the extremes of 
viscoelastic behavior to be expected from a 
"fiber-dominated" laminate layup versus a 
"matrix-dominated" laminate layup. 

• To perform creep tests using a multiple-angle 
laminate whose behavior was expected to be 
"fiber-dominated"; the layup selected was 
[ 0 / 45 / 0 /- 45 ] 6s - 

• To perform creep tests using a multiple-angle 
laminate whose behavior was expected to be 
"matrix-dominated"; the layup selected was 
[ 90 /- 45 / 90 / 45 ] 6s - 

• To develop a computer-based data reduction 
scheme which would provide whole-field strain 
contours based upon moire fringe patterns. 


TESTING EQUIPMENT 


Compression Fixture 

As previously mentioned, the compression fixture used in this 
study has been fully described in refs [12-14]. The fixture is a 
modification of the "end-loaded coupon" fixture, originally proposed 
by Clark and Lisagor! 15 !. The UW version of this device is shown 
schematically in Figurel. The fixture was designed to satisfy the 
following criteria: 

The fixture must: 

1 ) Apply compressive loading to the specimen using either: 

a) A SATEC lever arm creep frame, or 

b) A Tinius Olsen Universal testing frame. 

2) Insure that the applied load is "purely" uniaxial compression, 
applied through the centroid of the specimen cross-section. 

3) Rigidly constrain the specimen to prevent rigid body motion. 

4) Allow for easy inclusion in a moire interferometry optical 
system. 

Extensive fixture verification was conducted by Klein! 1 2 1 , through 
the use of strain gaged aluminum test specimens. It was found that 
the fixture does indeed produce a nearly uniaxial compressive load. 
The bending moments which are (inadvertently) induced result in 
bending strains whose magnitudes are typically less than 5% of the 
total axial strains.! 1 2 1 Hence, the effects of bending have been 
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ignored during the present study, i.e. it is assumed in all cases that 
the specimens are subjected to purely uniaxial compressive loading. 

Moire Fixture 

The moire interferometry fixture used in this research has been 
described in detail in ref [12-14]. The fixture uses a pair of two 
intersecting beams to produce a virtual reference grating. The 
virtual reference grating, in conjunction with a diffraction grating 
bonded to the specimen, results in a moir6 pattern. A sketch of the 
system is shown in Figure 2. A single 152 mm (6 inch) diameter 
laser beam is incident upon the fixture and mirrors A-D. These 
mirrors separate the incoming beam into four separate beams, which 
are then directed toward mirrors E-H. The second set of mirrors 
direct the beams of light to the specimen at an angle of 17.98°. This 
angle will produce a virtual reference grating with a frequency of 
1200 1/mm, when used with laser light produced by the argon-ion 
laser used in this study, i.e., light with a wavelength of 
514.4 nml 2 H 16 H 17 l. Mirrors E-H can be rotated in two out-of-plane 
directions, so as to precisely adjust the frequency and orientation of 
the virtual reference gratings. Mirrors A & C are used to create a 
virtual reference grating for use in producing moire patterns 
associated with "x H -displacements. Similarly, mirrors B & D are 
used to produce a virtual reference grating in the 

"y H -directionl 12 ]li8]. 
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Figure 3: Moire and Compression Fixture in the SATEC Creep Frame. 
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Figure 4: Moir6 Optical Set-up. 






12 


Testing Frames 

A SATEC Model G creep tester was used in all moire tests. 
This machine utilizes a dead weight lever arm system, with a 
maximum tensile load of 88,960 N (20,000 lbs). Figure 3 shows 
the moire and loading fixture mounted in the Satec creep frame. 

A Tinius Olsen universal load frame was used to determine the 
static ultimate compressive strength of the composite specimens. 
The Tinius Olsen can apply a maximum compressive load of 
266,880 N (60,000 lbs). 

>> 

Laser and Optics 

A Spectra Physics Series 2000 Argon Ion Laser, Model 
No. 2020-05, was used to produce a monochromatic light beam with 
a wavelength of 514.4 nm. The laser is capable of producing a laser 
beam with an intensity of up to 5 W, however only 200-300 mW 
were necessary to properly expose the photographic film used in this 
study. 

The light beam emanating from the laser head is expanded to 
form a 152 mm (6 in) dia beam of collimated light. The optical 
system used to expand and collimate the beam is shown 
schematically in Figure 4. The laser beam is first directed to a 
10X25 pm spatial filter using several beam steering mirrors. The 
light which passes through the spatial filter is diffracted (i.e. 
"bent"), causing the beam to expand in the cone-shaped envelope 
indicated in Figure 4. This expanding beam of light then impinges 
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upon a 152 mm dia parabolic mirror, which in turn reflects a 

collimated beam of light with the required 152 mm diameter. The 

parabolic mirror is adjusted so as to direct the reflected beam 

towards the moire fixture, providing the "incoming light beam" 
previously shown in Figure 2. 

Camera and Film 

A Cannon 35 mm camera, (base only), was used to photograph 
the moire fringe patterns, using shutter speeds ranging from 1/500 
to 1/2000 secs. Since the diffraction gratings used in this study 
had dimensions of 25.4 x 25.4 mm (1 x 1 in) the actual area of 

the moire pattern emanating from the specimen surface was 
645 mm 2 (1 in 2 ). A fl.4 lens was used to condense and focus the 
moire pattern on the 35 mm film plane. Kodak Tri-X Pan, black and 
white film was used to record the pattern. The film was processed 
according to the manufacturers recommended procedure. Prints of 
the moire patterns were produced using Kodak Poly Contrast Rapid II 
RC paper. Poly Contrast is a multi-contrast paper. A no. 4 contrast 
filter was used to produce satisfactory contrast. Standard 
development procedures were used to develop the prints. 

Specimen Preparation 

Nominal specimen dimensions were 152.4 mm (6 in) long, 
25.4 mm (1 in) wide and 6.35 mm (.25 in) thick. Approximately 
51 mm (2 in) of either end of the specimen is constrained within 
the end-loaded-coupon test fixture (see Figure 1), and hence the 
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unsupported gage length was nominally 51 mm (2 in). The Boeing 
Commercial Airplane Company supplied a IM6 graphite-epoxy [0]4g 
unidirectional panel and a IM6 graphite-epoxy [0/45/0/-45]6s panel. 
Imperial Chemical Industries (ICI) supplied an APC-2 graphite-PEEK 
[45/0/-45/90]6s panel. Specimens were cut from these three panels 
using a diamond coated abrasive disk. Holes (if required) were 
drilled using a carbide tipped drill. 

Moire grating molds, used to produce the specimen gratings, 
were prepared as described by Basehore and Post! 19 !. The specimen 
grating frequency was 600 l/mm. Crossed-gratings were used. 
Specimen preparation requires bonding the diffraction grating to the 
specimen surface. Specimen surface preparation is similar to 
techniques used when bonding resistance foil strain gages! 20 !. The 
moire grating mold is coated with a 500 A layer of aluminum using a 
vapor deposition system. Either Micro Measurement PC-10 or PC-6 
adhesive was used to bond the aluminum coating to the specimen 
surface. The adhesive is allowed to cure for 24 hours and excess 
epoxy is removed from the edges of the grating using a razor blade. 
The grating mold is then gently pried up, leaving the reflective 
diffraction grating bonded to the specimen surface. 

Test Matrix 

The selection of the laminate layup was conducted with the 
following in mind. 
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1) A determination of the fundamental engineering moduli Ei, E 2 , 
n -12 and 1121 was of interest. This information could be found 
by analysis of unidirectional specimens, both [0]48 and [90]48- 

2) Viscoelastic properties of multi-angle composite laminates 
were a concern. Miyano, Kanemitsu, Kunio and Kuhn! 21 ! have 
reported that some matrix dominated layups exhibit 
viscoelastic tendencies at room temperatures. Therefore, both 
a matrix dominated layup, [90/-45/90/45]6s. as well as a fiber 
dominated layup, [0/45/0/-45]6s. were investigated during 
this study. 

The [0]48 and [90]48 specimens were cut from a single 48-ply 
"parent" panel. The [0/45/0/-45]6s and [90/-45/90/45]6s specimens 
were cut from a single [0/45/0/-45]6s parent panel. 

Selection of Creep Stress Levels 

The quasi-static ultimate strength of each specimen was 
determined using specimens mounted in the Tinius-Olsen test frame. 
The ultimate strength tests are summarized in Table 1. The 
"specimen type" indicates both the material system and laminate 
layup (i.e., unidirectional or multi-angle laminate) for each 
specimen. The "specimen geometry" indicates the specimen hole 
diameter, if applicable. Each specimen was loaded in compression at 
a constant crosshead rate of 0.025 in/min. Audible cracking was 
heard throughout each test, up to and including specimen failure. 
The load at which complete catastrophic failure occurred was 
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recorded during each test and is listed as the "major failure load" in 
Table 1. The "major failure stress" was then calculated as the 
failure load divided by the net-cross sectional area, i.e., 
s fail = P/Anet- For specimens with holes a higher stress actually 
existed at failure, due to the stress-concentrating effect of the 
hole. No attempt was made to determine the actual maximum stress 
level induced in these specimens at failure. 

Once the ultimate strength of each specimen had been 
determined it was possible to select an appropriate creep stress 
level. The original intention was to test each specimen at a creep 
stress equal to 90% of the ultimate strength. However, in several 
tests it was found that a stress level of this magnitude resulted in 
moire fringe patterns which were so dense that individual fringes 
could not be identified. This was especially true of the 
"matrix-dominated" layups. Consequently creep stress levels were 
adjusted for each specimen so as to result in a moire pattern 
exhibiting a large number of fringes which could still be adequately 
interpreted. 

The creep stress levels used during the moire tests are 
summarized in Table 2. The "load applied" for each specimen type is 
indicated in the table. The corresponding "stress applied" equals the 
applied load divided by the net cross section. Again note that the 
concentrating effects of the holes are not accounted for in the value 
given for the applied stress. Finally, the "% of ultimate" relates the 


17 


creep stress level to the specimen ultimate strength, as previously 
presented in Table 1. 


TABLE 1 - ULTIMATE STRENGTH 


Test # 

Specimen Specimen 

type* Geometry 

Major Failure 
Load 

Major Failure 
Stress 

1 

A 

6.35 mm dia Hole 

77,800 N 

643 MPa 

2 

A 

No Hole 

102,300 N** 

634 MPa 

3 

B 

6.35 mm dia Hole 

64,500 N 

533 MPa 

4 

B 

No Hole 

102,300 N** 

634 MPa 

5 

C 

6.35 mm dia Hole 

24,500 N 

202 MPa 

6 

C 

No Hole 

33,400 N 

207 MPa 

7 

D 

6.35 mm dia Hole 

41,400 N 

342 MPa 

8 

D 

No Hole 

65,500 N 

400 MPa 

9 

E 

1.58 mm dia Hole 

65,600 N 

434 MPa 

1 0 

E 

3.17 mm dia Hole 

62,300 N 

441 MPa 

1 1 

E 

6.35 mm dia Hole 

40,000 N 

331 MPa 

1 2 

E 

9.52 mm dia Hole 

37,900 N 

375 MPa 


• Specimen types are defined as follows: 

Type A IM7-8551 graphite-epoxy [0]48 
Type B IM7-8551 graphite-epoxy [0/45/0/-45]6s 
Type C IM7-8551 graphite-epoxy [90]48 
Type D IM7-8551 graphite-epoxy [90/-45/90/45]6s 
TypeE APC-2 graphite-PEEK [45/0/-45/90] 6s 

The type A and B specimens with no hole had premature failures 
in the grip area. 


★ * 
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TABLE 2 - TEST MATRIX 

Specimen Specimen Load Stress 
Test # type* Geometry Applied Applied % of Ultimate 


1 

A 

6.35 

mm 

dia 

Hole 

40,000 

N 

330 

MPa 

62% 

2 

A 



No 

Hole 

80,100 

N 

496 

MPa 

78% 

3 

B 

6.35 

mm 

dia 

Hole 

50,700 

N 

419 

MPa 

78% 

4 

B 



No 

Hole 

66,700 

N 

413 

MPa 

65% 

5 

C 

6.35 

mm 

dia 

Hole 

8,000 

N 

66 

MPa 

33% 

6 

C 



No 

Hole 

5,300 

N 

33 

MPa 

16% 

7 

D 

6.35 

mm 

dia 

Hole 

18,700 

N 

154 

MPa 

45% 

8 

D 

6.35 

mm 

dia 

Hole 

37,400 

N 

308 

MPa 

90% 

9 

D 



No 

Hole 

18,700 

N 

115 

MPa 

29% 

1 0 

E 

1.58 

mm 

dia 

Hole 

50,700 

N 

335 

MPa 

77% 

1 1 

E 

3.17 

mm 

dia 

Hole 

42,700 

N 

302 

MPa 

69% 

1 2 

E 

6.35 

mm 

dia 

Hole 

37,400 

N 

308 

MPa 

93% 

1 3 

E 

9.52 

mm 

dia 

Hole 

24,000 

N 

238 

MPa 

63% 


* The specimen types are listed in Table 1. 


CALCULATION OF STRAIN FROM MOIRE FRINGES 

Moire fringe patterns represent lines of constant surface 
displacement. Although in some cases measured dispplacements are 
used directly! 22 ], normally the strain fields induced by external 
loading are of greater interest than the associated displacement 
fields. Hence, conversion of the displacements revealed by a moire 
fringe pattern to the corresponding strain field is usually desired. 
The moire patterns recorded during the present study were reduced 
to in-plane strain fields. The various techniques which have been 
used to perform this conversion can be roughly grouped into two 
categories; "mechanical differentiation"! 1 7 1 and the 
"displacement-field approach"! 23 !. The data reduction scheme used 
in the present study is based upon the displacement-field approach; 
both categories will be briefly explained below. 

Data Reduction by Mechanical Differentiation 

One category of reducing moire patterns to the corresponding 
strain fields includes methods based upon "mechanical 
differentiation" (also known as "moire of moire", "super moire", or 
"second-order moire" - the reader is warned that these terms are 
not uniformly applied, and different authors may use the same term 
to refer to a fundamentally different experimental technique). One 
way of performing mechanical differentiation is through the use of 
two identical transparencies (e.g., photographic negatives) of a 
moire fringe pattern. The two transparencies are initially exactly 
aligned, such that the original fringe pattern is observed. One of the 
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transparencies is then shifted with respect to the second by some 
amount, say Ax. This shifting results in a new "family" of fringes, 
which correspond to lines of constant (AN X /Ax). Hence, if the 
shifted Ax is relatively small, such that the quantity (AN X /Ax) 
approximates a true derivative, then these new fringes represent 
lines of constant strain, i.e., e x = 3u x /3x = (AN x /fAx), where 
f = virtual reference grating frequency A similar approach can be 
used to obtain the other derivatives required to characterize the 
in-plane strain field, specifically e y = 3v y /3y and 
©xy — (1/2)(du x /3y + 3v y /3x). 

Mechanical differentiation can be performed "directly" (i.e., 
using two transparencies as described in the preceding paragraph) or 
"indirectly". One indirect method is the use of double-exposure 
photography! 17 ]! 24 ]. In this case a single sheet of photographic film 
is mounted within a moveable film holder, and is exposed twice to 
the same moire fringe pattern. The film is shifted slightly between 
exposures. This approach is the optical equivalent to the physical 
shifting previously described, and results in a family of fringes 
representing constant strain. In a second indirect approach! 25 ! a 
digitized image of the original moire pattern is created and 
displayed on a CRT screen using appropriate computer-based 
hardware and software. A computer program then generates two 
versions of the digitized image, shifts one image with respect to the 
other, and displays both simultaneously. This procedure is the 
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numerical equivalent to the physical shifting of two transparencies, 
and again results in fringes of constant strain. 

An experimental methodology similar to mechanical 
differentiation is the use of a reference grating whose frequency 
has been adjusted so as to cancel a uniform portion of the 
displacement field. I 26 !. In this case the remaining fringes are 
related to local non-uniformities in the overall displacement field. 

Data Reduction by the Displacement-Field Approach 

Another category of reducing moire fringe patterns to strains 
includes methods based upon the "displacement-field approach"! 24 !. 
The technique is illustrated in Figure 5 and 6. A straight line which 
passes through a point of interest, such as line AB in Figure 5, is 
identified. The resulting displacements along this line are plotted. 
The slope of the displacement curve is the desired measure of 
strain, e.g. e x = 3u x /3x 3 (AN x /fAx). 

The data reduction scheme used in the present study is an 
automated version of the displacement-field approach illustrated in 
Figures 5 and 6. The technique was automated through the use of 
modern digitizing and computer equipment. The two primary pieces 
of equipment required were an AST Turboscan Digitizer and a 
Macintosh II computer with a 40 MB hard disk drive and a 40 MB tape 
backup. The moire patterns were first photographed with a 35 mm 
camera, and 203 x 254 mm (8x10 in) photographic prints were 
prepared, as described in a preceding section. Since the actual 
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Figure 5: Displacement-Position Plots for 9u/9x and 9u/9y.< 24 i 
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grating size was 25 x 25 mm (1 x 1 in), the photographic prints 
represent a nominal optical magnification of 8X. A digital record of 
the moire pattern is then obtained using the AST Turboscan 
Digitizer. This digitizer can be operated at a resolution ranging 
from 75 to 300 dots/in. In the present case the moire photos were 
scanned at a resolution of 300 dots/in, resulting in an effective 
resolution of about 2400 dots/in. The AST Turboscan Digitizing 
software allows for editing of the scanned image. The software 
permits viewing and editing of the image at the pixel level. This 
allows the user to "clean up" any imperfections (such as dust 
particles, scratches, etc.) which appear in the image. Once the 
image is cleaned to an acceptable level the image is stored on disk 
or tape backup. An example of the above procedure is shown in 
Figures 7 and 8. Figure 7 shows the original photograph of a moire 
pattern recorded for the y-displacements induced in an aluminum 
specimen with a hole, subjected to a 26 MPa (3770 PSI) 
compressive load. The specimen width and hole diameter are 
25.4 mm and 6.35 mm (1 in and 0.25 in) respectively. Although 
this is a high-quality moire image, note the various scratches and 
dust particles present. This image was scanned and "cleaned", 
resulting in the digital image shown in Figure 8. 
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Figure 7: V-Displacement Field, Aluminum Specimen. (Original) 
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Figure 8: V-Displacement Field, Aluminum Specimen. (Clean) 
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After the digitization and cleaning process was complete the 
moire data was reduced numerically using two computer programs 
developed in-house; programs "MOIRE" and "STRREG" (see 
Appendix A). Both programs were written using FORTRAN and 
standard Macintosh graphics software. Program MOIRE is used to 
convert the moir6 fringe patterns to strains. Program STRREG is 
then used to create a strain contour plot which can be displayed on a 
computer monitor or plotted using a graphics printer. Details of 
each program are given in the following two sections. 

Program MOIRE 

The digitized moire patterns created using the AST Turboscan 
system are stored in a bit map file, which in essence contains pixel 
information for every point on the image. The original 645 mm 2 
(1 in 2 ) grating area is represented by 5,760,000 pixel points in the 
digitized image. The program MOIRE computes derivatives in either 
the horizontal or vertical directions. Normal strain e x is 
determined from the u-displacement image, i.e. e x = du/dx, while 
e y is determined from the v-displacement image, i.e. e y = dv/dy. 
Note that the shear strain e xy can be found by taking the derivatives 
du/dy and dv/dx and summing them in the fashion: 
e xy = 1/2(du/dy + dv/dx). This latter calculation was not 
performed in this study. 

Program MOIRE performs the above calculations in the 
following order: 
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1 ) The current scan line (in the case of strain e x ) or scan column, 
(in the case of strain e y ) is identified. 

2) The data for the desired row or column is read from the image 
file. This data is stored as an 8 bit integer. The binary 
equivalent of this number represents a group of 8 pixels. A 
black pixel is stored with a 1 and a white pixel is stored with 
a 0. The 8 bit integer is converted to its binary equivalent and 
the pixel information is taken from each bit of this number and 
stored in a character array. 

2) The center of each fringe is located by scanning the array and 
locating the edges of each fringe. The center of the fringe is 
assumed to be midway between each of the two fringe edges. 

3) The appropriate derivative is calculated using the fringe 
center locations of three adjacent fringes (see Figure 9). For 
example, suppose the strain e x = du/dx is being calculated at 
fringe Nj, whose center is located at position Xj. The 
derivative is obtained using fringes immediately to the left 
and right of fringe Nj, i.e. fringes Nj_i and Nj+i, and so the 
derivative e x = du/dx can be approximated as: 

(n M - Nj.,)' 

f(*M ' X.-0. 




(Xi+ 1 , Ni+ 1 ) (Xi , Ni) (Xi- 1 , Ni- 1 ) 


Figure 9: Fringe Center Locations 
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where: 

f = virtual reference grating frequency 
Noting that: 

N i+1 - N M - ±2 
the above expression reduces to: 

I" 2* scale 

where: 

scale = scale factor added to account for the difference 
in size of the original specimen and the computer 
image. 

4) Strains are calculated at the center of each fringe along the 
entire row (or column), for every row (or column) in the image. 

5) The strain values, along with the corresponding x and y 
coordinates, are stored in a file for use during stage two. The 
maximum and minimum values determined throughout the 
image are also recorded in a file. 

Program MOIRE evaluates magnitudes of strain. The program 
does not use the actual fringe number when calculating strain, but 
rather relative fringe numbers that is, the fringe locations (Nj.-j and 
Nj+i) used to calculate strains e x or e y , at fringe number Nj are 
assigned values of 1 and 3 respectively. This process thus 
calculates the absolute value of the strain for the given strain field. 
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The strain fields in this study were either compressive (in the case 
of axial strains) or largely tensile (in the case of transverse 
strains). In order to analyze complex strain fields which contain 
both compressive and tensile strains, a more rigorous determination 
of local fringe order is required. 

Program STRREG 

The strain field map is generated using program STRREG. The 
strain map is generated using the file containing x and y 
coordinates and either e x or e y , and from the file containing the 
maximum and minimum strain values. 

1 ) The total range of strain, as determined using the maximum 
and minimum strain values, is broken up into eight intervals. 
A distinctive graphics pattern is assigned to each strain 
interval. 

2) Each calculated strain value and corresponding x- and 
y-coordinates are read from the strain file. 

3) STRREG opens a graphics window in memory, which is slightly 
larger than the original moire displacement photograph. The x- 
and y-coordinates are found in the graphics window. The 
corresponding graphics pattern is assigned to that location, 
and plotted to the next coordinate location. If the local strain 
value has increased or decreased into a different strain 
interval, the graphics pattern is changed accordingly. 

4) The strain field map is stored in a MacPaint™ format, allowing 
for easy editing and printing. 
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An example of the above approach is shown in Figure 10. The 
figure is based upon the original and digitized moire fringe patterns 
previously presented in Figures 7 and 8, respectively. 

Verification of the Data Reduction System 

The data reduction system was verified by comparing the 
results obtained for an aluminum specimen to the elasticity solution 
for an infinite plate with a circular hole! 24 !, and by comparison with 
a finite-element analysis for a finite-width plate. The elasticity 
solution is given by: 
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where 

s 0 = remotely applied uniform stress 
a = hole radius. 
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These equations were used to generate a data file containing x 
and y locations as well as axial strain values for an infinite plate 
with a hole. The STRREG program was then used to produce a strain 
field map. Figures 10 and 11 present the compressive axial strain 
fields determined for an aluminum plate using moire interferometry, 
the elasticity solution for an infinite plate, and a finite-element 
solution for a finite width plate, respectively. The comparison 
between measurement and theory is considered to be reasonable. The 
discrepancies which do exist, in the case of the elasticity solution, 
are due to the fact that this solution is valid for an infinite plate, 
while the experimental measurements were obtained using a 
specimen of finite size. The finite-element results, which assumes 
a finite width plate, shows a much better correlation. 
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strain 
E""5a strain 
strain 
B^B strain 


839 to 1533 ( ps/a ) 
1533 to 2226 (pa/a) 
2226 to 2920 (pa/a) 
2920 to 3613 (pa/a) 


strain • 3613 to *307 (pa/a) 
strain • *307 to 5000 (pa/a) 
strain ■ 5000 to 569* (pa/a) 
strain - 569* to 27890 (pa/a) 


Figure 10: Axial Strain Field Map for an Aluminum Specimen 
(Compressive Strains) 
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•train • 839 to 1533 (ps/a) 

l""5!51 strain - 1533 to 2226 (ps/a) 

mm strain « 2226 to 2920 (ps/a) 

mm strain • 2920 to 3613 (ps/a) 


stroin - 3613 to 430? (ps/s) 
atroin - 430? to 5000 (ps/s) 
strain • 5000 to 5694 (ps/s) 
strain - 5694 to 2?890 (ps/s) 


A: Infinite Plate (Elasticity Solution) 



rrm stroin - 839 to 1533 (ps/a) 

Stroin • 1533 to 2226 (ps/a) 

Wm strain - 2226 to 2920 (ps/a) 

■M strain - 2920 to 3613 (ps/a) 


stroin - 3613 to 4307 (po/s) 
stroin ■ 430? to 5000 (ps/s) 
strain ■ 5000 to 5694 (ps/a) 
strain ■ 5694 to 27890 (ps/s) 


B: Finite-Element 


Figure 1 1 : Axial Strain Field Maps Theoretical Solutions 
(Compressive Strains) 




RESULTS 


The results of the moire analysis will be presented in two 
separate sections below. Each test will be briefly discussed in an 
individual subsection. The tests involving specimens without holes 
will be described first, followed by a description of the tests 
involving specimens with holes. 

Specimens Without Holes 

Specimen type A: [0]48 Gr/Ep with no hole 

This specimen was tested at a creep stress level of -496 MPa 
(-72,000 psi), corresponding to 78% of ultimate stress. The strains 
measured for this specimen were used to calculate En and n 1 2 - 
Specifically, program MOIRE was used to calculate the axial and 
normal strains (ei and e 2 , respectively) induced throughout the gage 
section of the specimen. The average strain measurements were 
then determined to be: 

ei - -2982 pm/m 
e2 = 1632 pm/m 

These measurements correspond to the following values for the 
"major" Young's modulus and Poisson’s ratio: 

En = 166.3 GPa 
Hi 2 ■ 0.5473 

These values compare well with similar results previously reported 
for Gr/Ep material systems! 27 ]. Note the rather large value of Pi 2 . 
This high value is typical of Gr/Ep material systems. For these 
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materials Poisson's ratio is considerably higher in compression than 
in tension. For example, in Reference [28] values of n -12 ranging 
from about 0.3 to 0.25 are reported for tensile loadings, while 
values ranging from about 0.55 to 0.63 are reported for compressive 
loadings. 

Specimen type C: [90]48 Gr/Ep with no hole 

This specimen was initially tested at a creep stress level of - 
186 MPa (90% of the ultimate strength). This stress level resulted 
in a moire pattern which was far too dense for analysis. A reduction 
in stress level was therefore required; the stress was first reduced 
to -83 MPa and subsequently to -33 MPa (40% and 16% of ultimate 
strength, respectively). When the stress was reduced to this latter 
value a moire pattern with an acceptable fringe density was induced. 

The strains measured for this specimen were used to calculate 
E 22 and ri 2 i. The axial and normal strains (e 2 and e-i, respectively) 
induced throughout the gage section were determined using program 
MOIRE, and the average strain measurements were calculated as: 
e2 * -4143 pm/m 
ei - 138 pm/m 

These measurements correspond to the following values for the 
"minor" Young's modulus and Poisson's ratio: 

E 22 = 7.97 GPa 
ri2i = 0.0333 
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The measured value of E 22 corresponds well with that reported in 
Ref [28], 

The experimental measurements of En, E22, n-i2 and ri2i were 
evaluated to determine whether the following reciprocal 
relationship! 1 ] was satisfied: 
n2i = n 12 (E 2 2/Eii) 

Using the measured for En, E 22 and ni 2 listed above, the predicted 
value for the minor Poisson's ratio is ri 2 i » 0.0262. Hence, the 
measured value of 0.033 deviated from the "expected" value by about 
22%. This agreement is considered to be within the bounds of 
experimental error. The error represents the sum total of the errors 
in measuring each quantity, although it is likely due to a high 
percentage error in measuring ri2i, due to the very small transverse 
strains involved in the direct measurement of this quantity. 

The values of the measured moduli for the unidirectional 
graphite-epoxy specimens are summarized in Table 3. 

Specimen Type B: [0/45/0/-45]6 S Gr/Ep with no hole 

The specimen was tested at a creep stress level of -399 MPa 
(-57,900 psi), corresponding to 63% of ultimate strength. The axial 
and transverse strains (e y and e x , respectively) were determined 
using program MOIRE. The average values were: 
e y =* -4442 pm/m 
e x * 3769 pm/m 
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These measurements correspond to the following effective elastic 
moduli: 

Eyy = 89.8 GPa 
n yx = 0.8485 

These measured values were compared with the predictions obtained 
using a commercially-available computer program based upon 
classical lamination theory[29]. In the CLT analysis the measured 
values for En t E 22 and n 1 2 described above were assumed. The shear 
modulus G 12 was estimated as Gi 2 = 7.0GPa, based upon tabulated 
values in Ref [29]. The predicted effective moduli based upon the 
CLT analysis were: 

Eyy = 83.5 GPa 
Pyx = 0.9545 

Thus, the predicted and measured effective modulus and Poisson 
ratio differed by 7.5% and 11% respectively. 

Specimen Type D: [90/-45/90/45]6s Gr/Ep with no hole 

The specimen was tested at a creep stress level of -115 MPa 
(-16,700 psi), corresponding to 29% of ultimate strength. The axial 
and transverse strains (e y and e x , respectively) were determined 
using program MOIRE. The average values were: 
e y = -5382 jim/m 
e x = 906 pm/m 

These measurements correspond to the following effective elastic 
moduli: 

Eyy = 


21.4 GPa 
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n y x = 0.1683 

These measured values were compared with the predictions obtained 
using classical lamination theory. The predicted values were: 

E yy = 21.3 GPa 
n yx = 0.244 

Thus, the predicted effective modulus and Poisson ratio differed by 
0.5% and 31% from measured values, respectively. Note the rather 
large discrepancy between measured and predicted Poisson ratios in 
this case. The reason(s) for this discrepancy have not been 
identified. 

The comparisons between measured and predicted moduli for 
the multi-angle laminates are summarized in Table 3. 
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TABLE 3 

MEASURED AND PREDICTED MODULI 
FOR Gr/Ep LAMINATES 


Modulus 

Specimen type* 

Measured 

Moduli 

Predicted 

Moduli 

En 

Gr/Ep 

166.3 GPa 

NA 

E 22 

Gr/Ep 

7.97 GPa 

NA 

ni2 

Gr/Ep 

0.5473 

NA 

n2i 

Gr/Ep 

0.0333 

0.0262 

Eyy 

B 

89.8 GPa 

83.5 GPa 

X 

C 

B 

0.8485 

0.9545 

Eyy 

D 

21.4 GPa 

21.3 GPa 

riyx 

D 

0.1683 

0.244 

Specimen types 

are defined as 

follows: 



Gr/Ep IM7-8551 graphite-epoxy 

Type B IM7-8551 graphite-epoxy [0/45/0/-45]6s 

Type D IM7-8551 graphite-epoxy [90/-45/90/45]6s 


Specimens With Holes 

Specimen Type A: [0]48 Gr/Ep with 6.35 mm dia hole 

This specimen was tested at a creep stress level of -330 MPa 
(-47,900 psi), corresponding to 62% of ultimate strength. A 
"splitting" crack formed immediately upon application of the creep 
load. Although the crack was quite distinct and clearly visible, the 
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specimen supported the -330 MPa creep load, and hence was not 
considered to have "failed." The crack was very evident in the moire 
patterns, as shown in Figures 12 and 13. These photographs were 
taken 1 minute after application of the load. As indicated, the 
crack initiated at the top and bottom edges of the holes, and 
extended approximately 6.35 mm (.25 in) into the gage section in 
both directions. The creep stress was held constant throughout the 
10 hr test. No apparent crack extension occurred during the test; the 
moire patterns recorded after 10 hrs were essentially identical to 
those shown in Figures 12 and 13. 

The axial strain distributions calculated at a creep time of 1 
min and 10 hrs are presented in Figures 14 and 15 respectively. The 
crack had no effect on axial strain distribution, and is not at all 

apparent in these figures. No significant strain redistribution 
occurred during the 10 hr creep test. Axial strains were a maximum 
at the "3-o'clock" and "9-o'clock" positions around the hole, as would 
be expected. The overall distribution of axial strains is highly 
elongated in the fiber/loading direction. It is interesting to 

compare this axial strain distribution with the axial strain 

distribution for an isotropic aluminum plate, as presented in 

Figure 10. 

Equivalent results for the transverse strains are presented in 
Figures 16 and 17. In this case the existence of the crack is clearly 
evident in regions near the crack tip. Since the transverse strains 
were tensile, the crack tips were subjected to Mode I loading. 
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Observe that the transverse strain plots include a small "island" of 
high strain in the immediate vicinity of the crack tip. These islands 
are reminiscent of the isochromatic fringe "loops" normally seen for 
Mode I loading of a crack when using photoelastic techniques! 30 ]. No 
significant transverse strain redistribution occurred during the 10 
hr creep test. The u-displacement field (Figure 13) shows two 
regions (on the left and right sides of the hole) where the transverse 
strain changes sign, i.e., become compressive. As previously 

mentioned program MOIRE does not assign an algebraic sign to 
computed strains. The transverse strain field, therefore, does not 
indicate the small regions of compressive strain in the 

predominantly tensile strain field. 

A finite-element analysis was performed for this specimen 

geometry. The material properties used during the analysis have 
been previously listed in Table 3. The axial and transverse strain 
distributions predicted using the finite element analysis are 

presented in Figures 18 and 19, and should be compared with the 
experimental results in Figures 14 and 16 respectively. The 
comparison between measured and predicted axial strain 
distributions (Figures 14 and 18 respectively) is very good. Some 
differences in strain magnitudes are apparent, in particular near the 
regions of highest axial strains (the "3-o'clock" and "9-o'clock" 
positions). Nevertheless, the overall strain distributions agree very 
well. 
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Figure 12: V-displacement Field, [0]^, t = 1 min. 
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Figure 13: U-Displacement Field, [0] 48 , t = 1 min 


ORIGINAL PAGE lb 

©F POOR QUALITY 



Figure 14: Axial Strain Field Map [0]^, t * 1 min 
(Compressive Strains) 



strain • 471 to 845 (ps/s) strain 1 

strain - 845 to 1220 (pa/ a) — strain 

strain • 1220 to 1595 (pa/s) — strain 

strain • 1595 to 1969 (ps/s) BBB strain 


1969 to 2343 (ps/s) 

2343 to 2718 (ps/e) 

2718 to 3092 (ps/s) 

3092 to 24060 (ps/s) (Lagsnd is for both Figures) 


Figure 15: Axial Strain Field Map [0]^, t = 10 hrs 
(Compressive Strains) 
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Figure 16: Transverse Strain Field Map [0]^, t = 1 min 
(Tensile Strains) 



l: ! :': : :'-l strain - 450 to 785 (pa/s) — strain • 1792 to 2128 (pa/a) 

frrrrrri strain - 785 to 1121 (pa/s) — strain - 2128 to 2463 (pa/a) 

■M strain - 1121 to 1457 (pa/s) MB strain • 2463 to 2798 (po/a) t _ _ 

mm strain • 1457 to 1792 (pa/s) MB strain - 2798 to 13530 (pa/a) (Lagaod is for both Rguroa) 

Figure 17: Transverse Strain Field Map [0]^, t = 10 hrs 
(Tensile Strains) 
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L...I strain - 171 to 815 (po/o) 

t^.wJ strain - SIS to 1220 (po/o) 

■■■ strain • 1220 to 1595 (po/o) 

■■■ strain - 1595 to 1969 (po/o) 


strain • 1969 to 2313 (po/o) 
strain • 2313 to 2718 (po/o) 
strain • 2718 to 3092 (po/o) 
strain • 3092 to 21060 (po/o) 


Figure 18: F.E. Axial Strain Field Map [0]^ 
(Compressive Strains +) 



ES3 strain - 70S to 1121 (po/o) M strain - 2128 to 2163 (po/o) 

— strain • 1121 to 1157 (po/o) strain - 2163 to 2798 (po/o) 

strain • 1157 to 1792 (po/o) ■■■ strain • 2798 to 13530 (po/o) 


Figure 19: F.E. Transverse Strain Field Map [0]^ 
(Tensile Strains +) 
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On the other hand, the comparison between measured and 
predicted transverse strains (Figures 16 and 19, respectively) is 
not so good. The difference is due to the existence of the crack. The 
predicted distribution shows a region of high transverse strain in 
the "12-o’clock" and "6-o'clock" positions. This is precisely where 
the splitting crack originated during the test. The existence of the 
crack unloaded material in the immediate vicinity of the hole, 
reducing transverse strains near the hole. Strains were 
redistributed such that a general increase in transverse strain at 
regions removed from the hole occurred. The finite element analysis 
predicts several regions of compressive strain, on the right and left 
sides of the hole, and also at the edge of the hole on ±45° diagonals. 
The experimental u-displacement field shows the two regions on the 
left and right sides of the hole, but the regions on the ±45° diagonals 
are not apparent. 

Specimen Type C: [90]48 Gr/Ep with 6.35 mm dia hole 

This specimen was tested at a creep stress level of -66 MPa 
(33% of the ultimate strength). During the test severe grating 
degradation took place, resulting in excessive "noise" in the 10 hour 
strain field maps. The degradation resulted in an overall loss of 
reflection and refraction of the grating; the specimen grating 
appears to "cloud over". It is suspected that the degradation is due 
to oxidation of the grating. The phenomena occurs randomly, and 
was also occasionally encountered by Klein! 1 2]. in the present test 
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the overall shape of the strain distributions was still evident, 
despite the degraded moire pattern. 

The axial strain distributions calculated at a creep time of 
1 min and lOhrs are presented in Figures 20 and 21 respectively. 
An analysis of the two figures shows that some strain 
redistribution occurred during the 10 hr creep test. This 
redistribution is most clearly evident in regions near the 
"12-o'clock" and "6-o'clock" positions. At time t = 1 min relatively 
large regions of small axial strain are evident at these locations, 
while at t= lOhrs these regions have been significantly reduced 
in size. Upon further inspection the broad "winged" shaped regions, 
roughly defined by the ±45° directions, are seen to have generally 
expanded during the 10-hr creep test, indicating axial strain 
redistribution due to viscoelastic effects. Maximum axial strains 
occur at the sides of the hole, as in the type A specimen. 

The transverse strain fields for the 1 min and 10 hr creep 
times are presented in Figures 22 and 23. The transverse strain 

fields also show strain redistribution during the 10 hr creep test. 
Notice that "loops" of high strain exist at the "3-, 6-, 9- and 12- 
o'clock" position around the hole. In the 10 hr strain field, these 
loops of higher strain have been expanded, indicating viscoelastic 
behavior. 

The axial and transverse strain distributions predicted using 
the finite-element analysis are presented in Figures 24 and 25, and 
should be compared with the experimental results in Figures 20 and 
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22, respectively. The overall comparison between predicted and 
measured axial strains is considered to be reasonable, although 
differences in strain magnitudes are apparent. The finite-element 
prediction for the transverse strain field however, is very different 
from the strain field measured experimentally. The primary 
characteristic of the finite-element prediction is the highly 
elongated "lobes" of transverse strain, centered symmetrically 
above and below the hole. Also note that the predicted transverse 
strains are roughly one order of magnitude less than the measured 
transverse strains. In addition, the finite-element transverse strain 
field shows small regions of compressive strain at the edge of the 
hole on ±45° diagonals, as observed for the [0]48 specimen. These 
small regions of compressive strain were not apparent in the 
experimentally-determined transverse strain field. 

The discrepancies between prediction and measurement 
evidenced in this case is far greater than the discrepancies 
associated with any of the other comparisons considered during this 
study. The reason(s) for this isolated, but severe, discrepancy have 
not been identified. Perhaps the test specimen developed some 
significant internal defect upon application of the creep load, such 
as a large internal delamination, for example. Such a hypothesized 
defect was not apparent, however, and it is likely that if such a 
defect existed it would have been readily evident in the moire 
pattern. 
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strain • 5913 to 7053 (pa/a) 

strain ■ 7053 to 9193 (pa/a) 

strain • 9193 to 9333 (pa/a) , . . 

strain - 9333 to 19450 (pa/a) (Logand is for both Figures) 


<: • strain • 1352 to 2192 (pa/s) 

ESE3 strain ■ 2492 to 3632 (pa/s) 

■■i strain - 3632 to 4772 (pa/s) 

■IH strain « 4772 to 5913 (pa/s) 


Figure 21 : Axial Strain Field Map [90]^, t = 10 hrs 
(Compressive Strains) 
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Figure 22: Transverse Strain Field Map [90] 48 , t = 1 min 
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Figure 23: Transverse Strain Field Map [90]^, t = 10 hrs 
(Tensile Strains) 
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■■1 strain • 3632 to 4772 (pa/a) 

■■I strain • 4772 to 5913 (pa/a) 


strain • 5913 to 7053 (pa/a) 
strain ■ 7053 to 6193 (pa/a) 
stroln ■ 0193 to 9333 (pa/a) 
strain - 9333 to 10450 (po/a) 


Figure 24: F.E. Axial Strain Field Map [90]^ 
(Compressive Strains +) 
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Figure 25: F.E. Transverse Strain Field Map [90]^ 
(Tensile Strains +) 
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Specimen Type B: [0/45/0/-45]6s Gr/Ep with 6.35 mm dia 
hole 

This specimen was tested at a creep stress level of -419 MPa 
(-60,770 psi), corresponding to 78% of ultimate strength. The 

moire patterns associated with the v- and u-displacements obtained 
at t = 1 min are shown in Figures 26 and 27. Nearly identical u- 
and v-displacement fields for quasi-isotropic specimens with holes 
under tensile loading were reported by Postl 2 !. 

The axial strain field maps, shown in Figures 28 and 29, show 

no significant redistribution of strain during the 10 hr creep test. It 

is interesting to note that the disturbances in the strain field 
resulting from the hole, in the axial strain field do not extend as far 
from the hole as in the [0]48 (type A) axial strain fields (Figures 14 
and 15). The type B axial strain fields are similar in shape to those 
of the isotropic aluminum plate, shown in Figure 10. 

The transverse strain fields for this specimen are presented in 
Figures 30 and 31. These strain fields also show no signs of 

viscoelastic behavior for the 10 hr creep test. Note that areas of 
high strain exist at the top and bottom of the hole as well as along 
±45° diagonals extending away from the hole. Klein found similar 
double strain "peaks" in his analysis [Ref 12, Figure 28]. 

Finite-element strain field predictions are presented in 
Figures 32 and 33. A comparison of these results with the 

experimental results shows a very good overall correlation in both 
the axial and transverse strain distributions. However, the 
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finite-element transverse strain field shows small regions of 
compressive strain, not observed in the experimental strain field, at 
the edge of the hole on ±45° diagonals, just as in the [0]48 specimen. 

Specimen Type D: [90/-45/90/45]6s Gr/Ep with 6.35 mm dia 
hole 

The specimen was initially tested at a creep stress level of 
-308 MPa (-44,671 psi), corresponding to 90% of the ultimate 
strength. This stress level resulted in an immediate delamination of 
the specimen on both sides of the hole. This delamination slowly 
grew with time. Figures 34 and 35 show the delamination growth 
for the specimen. Strain field maps were not generated for this 
specimen due to the extremely high fringe densities resulting from 
the delamination. Complete, catastrophic failure of the specimen 
occurred at about 16 hours after loading. The failed specimen is 
shown in Figures 36 and 37. The failed region of this specimen 
resembles the brooming failure observed for the ultimate strength 
test specimens. Note the relatively symmetric failure mode, 
indicating an absence of overall buckling failure. The growth of the 
delamination throughout the 16 hr test indicates viscoelastic 
behavior for this layup at room temperature, and also clearly shows 
the intimate relationship between viscoelastic and fracture 
mechanisms. Similar failures of a composite specimen subjected to 
fatigue loadings have been reported by Carlssonl 32 !. 

A second type D specimen was prepared and this specimen was 
tested at a creep stress level of -154 MPa (-22,335 psi), 45% of 
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the ultimate strength. The purpose of the second test was to see if 
viscoelastic effects would be evident at this lower stress level. 

Figures 38 and 39 present the axial strain fields. There are 
signs of slight strain redistribution for this test; the level of 
redistribution is on the same order as that presented in Figures 20 
through 23 for the type C [90]48 specimen. It is important to note 
that viscoelastic behavior took place at room temperature over a 
10 hr test period, at a stress level that was less than half of the 
ultimate stress. The shape of the strain fields are similar to those 
presented in the type B specimen layup. 

Slight viscoelastic response can also be seen when comparing 
the 1 min and 10 hr transverse strain fields, presented in 
Figures 40 and 41. The transverse strain field shows the same 
general shape as the type B specimen with similar areas of strain 
concentrations. The transverse strain fields show areas of "jagged" 
edges, this results from the relatively few fringes in the 
u-displacement field and hence few data points. Comparisons of 
the finite-element strain field predictions, shown in Figures 42 and 
43, to the experimental strain fields show very similar strain 
distributions in both the axial and transverse cases. However, the 
finite-element predicts small regions of compressive strain, at the 
edge of the hole on ±45° diagonals, just as in the other specimens. 
These areas of compressive strain were not apparent in the moire 
patterns. 
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Figure 29: Axial Strain Field Map [0/45/0/-45) 6s , t = 10 hrs 
(Compressive Strains) 


Figure 28: Axial Strain Field Map [0/45/0/-45],,, t = 1 min 
(Compressive Strains) 
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Figure 30: Transverse Strain Field Map [0/45/07-45]^, t = 1 min 

(Tensile Strains) 
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strain - 1779 to 18150 (po/o) (Loflond is for both Figures) 


Figure 31 : Transverse Strain Field Map [0/45/07-45]^, t = 10 hrs 

(Tensile Strains) 
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Figure 32: F.E. Axial Strain Field Map [0/45/0/- 45]^ 
(Compressive Strains +) 


r 1 ' ' i 



0 (pa/a) 
1360 (pa/a) 
1930 (pa/a) 
2499 (pa/a) 
3069 (pa/o) 


stroin - 3069 to 3639 (pa/a) 
strain • 3639 to 4209 (pa/a) 
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Figure 33: F.E. Transverse Strain Field Map [OMS/OMS]^ 
(Tensile Strains +) 
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TIME = 1 HRS TIME = 5 HRS TIME = 1 0 HRS 


Figure 34: Delamination of Type D [90/45/90/-45] 6s Gr/Ep V-Displacement 
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Figure 35: Delamination of Type D [90/45/90/-45] 6s Gr/Ep U-displacement 
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Figure 37: Side View of Failed [90%45790°/45°] 6s Gr/Ep Specimen. 



ORIGINai 


Rage IS 


0F P °OR QUALITY 



Figure 38: Axial Strain Field Map [90/-45/90/45],,,, t = 1 min 
(Compressive Strains) 
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Figure 39: Axial Strain Field Map [90/-45/90/45] 6s , t = 10 hrs 
(Compressive Strains) 




Figure 40: Transverse Strain Field Map [90/-45/90/45] 68 , t = 

(Tensile Strains) 
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Figure 41 : Transverse Strain Field Map [90/-45/90/45] 64 , t 

(Tensile Strains) 
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(Legend is for both Figures) 
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Figure 42: F.E. Axial Strain Field Map [907-45/90/45]^ 
(Compressive Strains +) 
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Figure 43: F.E. Transverse Strain Field Map [90/-45/90/45], 

(Tensile Strains +) 
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Graphite/PEEK Specimens with holes 

The Gr/PEEK specimens tested were quasi-isotropic Type E: 
[45/0/-45/90]6s- The four specimens tested contained holes of the 
following diameters; 1.58 mm, 3.17 mm, 6.35 mm and 9.52 mm. 
The stress levels of the test were: -335 MPa (-48,588 psi), 

corresponding to 77% of the ultimate strength; -302 MPa 
(-43,800 psi), corresponding to 69% of the ultimate strength; 
-308 MPa (-44,671 psi), corresponding to 93% of the ultimate 
strength; and -238 MPa (-34,518 psi), corresponding to 63% of the 
ultimate strength; respectively. No viscoelastic response was found 
in either the axial or transverse strain distributions for any of the 
four Gr/PEEK specimens tested, and therefore only the 1 min strain 
fields will be presented. 

The axial strain field maps are presented in Figures 44 through 
47. The axial strain field for the specimen with a 1.58 mm dia hole 
(Figure 44) shows a near uniform strain field in a significant 
portion of the test area. The axial strain field for the 3.17 mm dia 
hole specimen (Figure 45) shows the familiar "wing" shaped strain 
regions found in the other specimens. The disturbances generated by 
the hole extend further from the center of the strain field. The 
6.35 mm hole specimen shows an axial strain field (Figure 46) 
similar to those of the 6.35 mm hole Gr/Ep specimens. The 
disruptions in the strain field, in the familiar "winged" shaped 
regions, are much more extensive than in the specimens with 
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smaller holes. The disturbances extend throughout the entire gage 
area of the specimen. 

The axial strain field map for the specimen with a 9.52 dia 
hole is presented in Figure 47. The strain field shows some of the 
"wing" shaped strain regions, however, most of the strain patterns 
extend out past the gage area to the edge of the specimen where they 
combine with effects produced by the edge of the specimen. 

The transverse strain fields are presented in Figures 48 
through 51. The 1.58 mm dia hole specimen transverse strain field 
(Figure 48) shows large portions of uniform strain. The outer edges 
of the gage area show almost no effects from the hole. The general 
shape of the strain field near the hole however, is similar to the 

strain fields found in the other specimens. Small "looping" regions 

of strain emanate on +45° diagonals from the hole. The transverse 
strain fields for the 3.17 mm dia hole specimen (Figure 49) shows 
an increase in these "looping" regions. The diagonal direction loops 
in the 6.35 mm dia hole specimen (Figure 50) have extended to the 
corners of the gage area. It is also interesting to note how the edge 
effects have further combined with effects from the hole. The 

general shape of the transverse strain, like the axial strain field for 
this specimen, is similar to the shapes found in the Gr/Ep 

specimens. The transverse strain field map for the 9.52 mm dia 
hole specimen is presented in Figure 51. The "looping" strain 
regions are very apparent in this specimen. Several areas of the 
strain field display regions where a strain pattern extends away 
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Figure 45: Axial Strain Field Map, 1/8" dia hole Gr/PEEK, t = 1 min 

(Compressive Strains) 
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Figure 44: Axial Strain Field Map, 1/16" dia hole Gr/PEEK, t = 

(Compressive Strains) 
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Figure 46: Axial Strain Field Map, 1/4" dia hole Gr/PEEK, t = 

(Compressive Strains) 
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Figure 47: Axial Strain Field Map, 3/8" dia hole Gr/PEEK, t = 1 min 

(Compressive Strains) 
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Figure 48: Transverse Strain Field Map, 1/16" dia hole Gr/PEEK, t = 1 min 

(Tensile Strains) 
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Figure 49: Transverse Strain Field Map, 1/8" dia hole Gr/PEEK, t = 1 min 

(Tensile Strains) 
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Figure 50: Transverse Strain Field Map, 1/4" dia hole Gr/PEEK, t = 1 min 
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Figure 51 : Transverse Strain Field Map, 3/8" dia hole Gr/PEEK, t = 1 min 

(Tensile Strains) 
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Figure 52: F.E. Axial Strain Field Map, 1/4" dia hole Gr/PEEK 
(Compressive Strains +) 
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Figure 53: F.E. Transverse Strain Field Map, 1/4" dia hole Gr/PEEK 

(Tensile Strains +) 
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from the hole and then "loops" back on itself, ending where it 
originated. 

A finite-element analysis was performed for the 6.35 mm dia 
hole specimen using the following material properties reported by 
Tsai[30] : 

En = 134 GPa 
E 22 = 8.9 GPa 

G 12 = 5.1 GPa 

A value of n-i 2 = 0.28 was reported by Tsai for Gr/PEEK in tension. 
Since n 1 2 has been found to be significantly higher in compression, 
than in tension for other composite systems, a value of n-12 = 0.5 
was arbitrarily selected for use in the finite-element analysis. 
Comparison of the finite-element strain fields (Figures 52 and 53) 
with the experimentally measured strain fields (Figures 44 and 48) 
shows a good correlation between the two. 


SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


Summary 

The overall goal of this study was to measure elastic and 
viscoelastic strain fields induced in composite laminates by 
compressive loading. Moire interferometry was used to produce 
moire fringe patterns for the u- and v-displacement fields. A 
computer based data reduction system was developed which, using a 
digitized image of the u- and v-displacement field, produced a whole 
field strain field map for both the transverse and axial strain fields. 
It would also be possible to use a similar approach to produce a 
shear strain field map, although this was not attempted in this 
study. These strain field maps allowed for comparison of strain 
distributions at different times. A comparison was made for strain 
fields of different composite material layups and also for strain 
fields in composites specimens with different size holes. Solid 
specimens were tested and the average strain computed using the 
computer data reduction system. These strains were then used to 
determine fundamental engineering moduli. Measured strain 
distributions for specimens with holes were compared with 
predictions obtained using finite-element techniques. 
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Conclusions 
Strain Field Map 

The strain field map proved to be a valuable tool in the 
analysis of strain distributions. Comparisons of different material 
layups showed distinctly different stain fields. The irregularities 
that appeared in the strain fields did not mask the overall shape of 
the strain fields, so analysis and comparisons could still be made. 
The "irregularities" in the experimentally measured strain 
distributions might at first be interpreted as "noise" or "scatter" in 
the data. While a percentage of these irregularities is undoubtedly 
due to experimental error, the authors are convinced that a 
significant fraction of this scatter is, in fact, an accurate 
representation of the surface strains. That is, the "scatter" in 
strain is simply a reflection of the inherent heterogeneity of 
composites. This is supported by the results of Post and others! 6 !, 
where composite specimens subjected to tensile loadings have 
exhibited irregularities in surface strains as high as ±15%. The 
nonuniformity in surface strain distributions tend to cause 
irregularities in the strain contour maps presented herein. Note that 
the finite-element analysis does not account for the possibility of 
inhomogeneities, and hence the strain contours prepared using 
finite-element results indicate very smooth and well-defined strain 
contours. 

It is important to note once again that the data reduction 
scheme used in this analysis measures magnitudes of strain. The 
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[ 0 ] 4 8 transverse displacement field showed two "islands" of 
compressive strain in a predominantly tensile field. The data 
reduction method does not identify the algebraic sign of strain. The 
responsibility falls on the user to identify any regions where sign 
changes may occur. The strain fields analyzed in this study were 
relatively simple, so this could be done easily. However, if complex 
strain fields are to be analyzed using this approach, a more rigorous 
determination of actual fringe number is required. 

Viscoelastic Effects 

Viscoelastic effects were noted in two of the specimens 
containing holes. The matrix dominated layups, [90]48 and 
[90/-45/90/45]6s» showed mild viscoelastic response over the 
length of the test when subjected to stress levels of about 33% and 
45% of the static ultimate strength, respectively. Major time 
dependent effects were seen in a [90/-45/90/45j6s specimen with 
hole subjected to a creep load of -308 MPa, or about 90% of the 
static ultimate strength. This specimen delaminated immediately 
upon application of the load, and the delaminated region grew until 
complete catastrophic failure occurred at about 16 hours after 
application of the load. In this case the response is a complex 
combination of both viscoelastic and fracture mechanisms. 

With the exception of the [90/-45/90/45]6 S specimen which 
failed catastrophically after 16 hrs, the viscoelastic effects 
observed during this study were relatively mild. Note, however, that 
the measurements were obtained over a 10-hr period and at room 
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temperatures. The implication is that at much longer times (say 
after weeks, month, or years) the composite laminates considered in 
the study may exhibit very significant viscoelastic deformation, 
which may be accompanied by delaminations and eventually 
catastrophic failure. 

It is believed that the (nominal) 10-hr creep test time used in 
this study is "about" as long in duration as can be conveniently 
performed using moire interferometry. Hence, if truely long-term 
effects are to be investigated using moire interferometry, the 
effects due to viscoelastic behavior must be "accelerated" in some 
fashion. A recommendation for further study is use of the present 
moire system under conditions of moderately elevated temperatures 
(say 90°C, or about 200°F). 

Comparison Between Experimental and Finite-Element 

Results 

The finite-element resuits compare well in most cases with 
the experimental results. However, in general the finite-element 

strain field maps indicate lower strain in areas removed from the 
hole and higher strain in areas close to the hole. It is hypothesized 
that the high stress levels at the edge of the hole induced nonlinear 
material behavior at these points. This nonlinear behavior may arise 
from the formation of localized micro-cracks or voids. These micro 
failures in the vicinity of the hole would result in a redistribution 
of strain. This redistribution may result in slightly lower strains 
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immediately adjacent to the hole and slightly higher strains 
throughout the rest of the strain field, which agrees in a qualitative 
sense with the results obtained. 

Overall Conclusions regarding the Use of Moire 

Interferometry 

To the author's knowledge, this program represents the first 
time that moire interferometry has been used under creep 
conditions. Overall, the use of moire interferometry has been very 
successful. The moir6 system developed has proven to be stable and 
suitable for use during creep periods up to 24 hrs in duration. It is 
not clear if the system would remain stable for much longer times, 
say for several days or weeks. It has been possible to identify 
slight viscoelastic deformations which occur at room temperatures 
in composite laminates through the use of moire. It is not likely that 
other popular strain measurement devices such as resistance foil 
strain gages would have been stable enough to have identified this 
behavior. 

Moire interferometry itself is a powerful tool for use in studying 
composites. An outstanding advantage of moire is the ability to 
obtain whole field displacements and/or strains, as opposed to 
point-measurement techniques such as strain gages, extensometers, 
or LVDT's. Furthermore, moire interferometry can be configured so 
as to obtain measurements with a sensitivity which equals or 
exceeds those available with other techniques. 
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On the other hand, moire interferometry is an exacting 
technique. Specimen and grating preparation techniques, as well as 
subsequent bonding of the grating to the specimen, are considerably 
more involved than those associated with bonding strain gages, for 
example. The personnel involved in this study have found that as 
more experience was gained, much better moire patterns were 

obtained. These improvements occurred even though no apparent 

changes in procedure were implemented. In other words, 
considerable expertise is required before high-quality moire 
patterns can be obtained by the researcher. This expertise can only 
be gained through personal experience in the lab. ‘ The exacting 
nature of moire interferometry has definite practical implications. 

For example, in most experimental studies involving the use of 

strain gages several repeated tests, over a wide range of test 
conditions, will normally be performed. Very few repeated tests 
were performed during this study, and the range of test conditions 
was more limited than might otherwise be desired. These aspects of 
the program were required by the combination of time and budget 
constraints and the degree of difficulty in performing each moire 
test. 

Recommendations 
Moire Fixture 

The Moire fixture could be improved with additional alignment 
features, which would insure the specimen face was perpendicular 
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to the incoming light and centered in the loading frame. The ability 
to use multiple frequencies of virtual gratings would also aid in the 
examination of high strain fields. In applications such as the matrix 
dominated layups, a different frequency of grating for the 
transverse and axial strains would greatly improve the analysis of 
these strain fields. This may be achieved using the 
recently-developed "composite" specimen gratings described by 
Wang, Kang and KobayashiP 2 ). 


Elevated Temperature/Humidities Testing 

Viscoelastic response was seen in some of the specimens 
tested. It would be of great interest to expose the specimens to a 
more aggressive environment such as elevated temperatures and 
humidities where viscoelastic effects are likely to over a 
convenient test period. This type of testing would allow for 

viscoelastic analysis of any material layup such as the 
"fiber-dominated" layups which do not show any viscoelastic 
response at short times and at ambient conditions. 

Elevated temperature/humidity testing would require 
additional apparatus in the form of an environmental chamber to fit 
around the compression and/or moire fixture. However, the 
compression fixture and moire system developed during this study 
were designed with eventual inclusion of an environmental chamber 
in mind. Hence, the inclusion of a chamber in the existing test setup 
can be achieved with relative ease. 
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Data Reduction System 

The data reduction system developed during this study resulted 
in generally good results. However when moire fringe patterns of 
low quality were obtained, i.e., excessive noise was present, the 
resulting strain fields also showed areas of noise. Solutions for 
this are of two types. Either obtain consistently "clean” moire 
patterns, or develop a computer-based filtering system which would 
allow for a numerical filtering of (digitized) "noisy" moire patterns 
to obtain clean patterns. Potential methods of numerical filtering 
of the digitized moire patterns were briefly evaluated during this 
study, but a satisfactory system has not yet been developed. 
Numerical filtering should be possible, however, and is recommended 
as an area of further study. 

As previously mentioned, as experience is gained by the 
individual the quality of moire patterns obtained by that person 
generally improves. Hence, the only way of obtaining consistently 
"good" moire patterns is the development of personal experience. 

A second area of improvement is in calculation of strain 
values. In the present effort strains were determined using a simple 
averaging scheme which approximates the derivative of a 
displacement field. An improvement would be to incorporate a more 
sophisticated numerical differentiation scheme, such as the use of 
Lagragian polynomials. 
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APPENDIX A:COMPUTER PROGRAMS 


Two FORTRAN-based computer programs were developed 
in-house during this program; programs MOIRE and STRREG. These 
programs allowed the calculation of strains from digitized moire 
fringe patterns, and subsequent display of strain contour plots. 
Listing of the source code for each of these programs is listed 
below. 

Program MOIRE 


implicit none 
integer*2 data(6) 
integer*2 width .height 
integers pic(400) 
integer*1 testbyte.testnum 
real df, scale, evencheck 
real str(IOOO) 
character‘1 pix(8) 
integer Isf, rsf ,tof ,bof , coif, coll 
real 

f ri ng e( 1 000), ddymin.ddy max, ddxmin,ddxmax,tstr,avestr, rent 
integer resp.cnt 
character filename*64 
integer*1 offset 

integer i, k, j, n, I, m, scancolm, scan line, row, col, fileloc, byte 
integer sec 

character^ pixel(8,3000) 

tstr = 0.0 
ent = 0 


90 

ddxmin = 1.0 
ddxmax = -1.0 


ddymin = 1 .0 
ddymax = -1.0 


* Open the file. Bitmap data is binary, so use unformatted access. 

write(9,*)'lnput the file to be scanned' 
read(9,500)filename 

500 format(A) 

open(unit=20, FILE=filename, STATUS= , old , ,ACCESS='direct', 
+ FORM='unformatted',RECL=2) 

write(9,501 ) 

501 format(1x, 'Choose the desired analysis', 

+ /1x,'1. Calculation of du/dx', 

+ /lx, '2. Calculation of dv/dy',/1x) 
read(9,*)resp 

write(9,*)'lnput the scale. (Pixels/In).' 
read(9,*)scale 

write(9,*)" 

write(9,*)'lnput the reference grating frequency. (I/mm)' 
read(9,*)df 

df = 1.0/(df*25.4) 

write(9,*)'reading the moire picture' 

★ 

* read file info 
do (i = 1 ,6) 

read(unit=20,REC=i)data(i) 

repeat 

close(20) 
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width = data(4) 
height = data(5) 

write(9,*)’scanning resolution \data(1) 
write(9,*)'scanning width \8*data(4) 

write(9,*)'scanning height \data(5) 
write(9,*)' ' 

* 

* This block checks to see if the scanning width is even or odd 

* If the width is even an extra byte is added to each record in the 

* scanned file, so this must be compenstated for here. 

* 


evencheck = width/2.0 
evencheck = evencheck - (width/2) 
if(evencheck.EQ..5)then 
offset = 3 

else 

offset = 4 

endif 


* This block inserts the data into a screen size array 


IF(resp.EQ.1)THEN 

open(unit«10,file= , ddx.output' > status= , new') 
do(scanline = 40,height,8) 
sec = (scanline - 1 + 8)/8 
write(9,*)'Scanning Section # \sec 

open(unit=20, FILE=filename, STATUS='old\ACCESS='direct\ 
+ FORM=’unformatted\RECL=1) 
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do(row = 4,8,4) 

fileloc=(row-1 +scanline-1 )*(width+offset)+51 3 

do(byte=fileloc,fileloc+width+offset) 

i=byte-fileloc+1 

read(unit=20,rec=byte,end=200)pic(i) 

repeat 

200 do(i=1 , width) 
k = i + 4 

if(pic(k).LT.O)then 
pix(1) = 'B' 

else 

pix(1) = 'W 

endif 

scancolm = ((i*8)-7) 
pixel(row,scancolm) = pix(1) 
testnum = 64; 

do(j=2,8) 

testbyte = pic(k). and. testnum 
if(testbyte.NE.O)then 
pix(j) » -B 1 

else 

pix(j) = W 

endif 

testnum = testnum/2 

scancolm = ((i*8)-8+j) 

pixel(row,scancolm)=pix(j) 

repeat 

repeat 

repeat 


close(20) 
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CALL 

D DX(width, height, pixel, scale, df.ddxmin.ddxmax, sec, tstr.cnt) 
repeat 

ELSEIF(resp.EQ.2)then 


open(unit=1 2, file='ddy. output', status='new’) 
do(sec * 1, width) 

write(9,*)'Scanning Section # \sec 

open(unit=20, FILE=filename, STATUS= , old' I ACCESS='direct’, 
+ FORM»'unformatted',RECL=1) 

do(row = 40, height) 

ftleloc=(51 2+offset) + (sec-1) + (row-1 )*(width+offset) 
byte=fileloc 

read(unit=20,rec=byte,end=21 0)pic(1 ) 

210 if(pic(1).LT.0)then 
pix(1) = 'B' 

else 

pix(1) = ’W’ 

endif 

scancolm = (1) 
pixel(scancolm.row) = pix(1 ) 
testnum = 64; 

do(j=2,8) 

testbyte = pic(1).and.testnum 
if(testbyte.NE.O)then 
pix(j) = 'B' 

else 


pix(j) = 'W' 
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endif 

testnum = testnum/2 

scancolm = (j) 

pixel(scancolm,row)=pix(j) 

repeat 

repeat 

close(20) 

CALL 

D D Y(width, height, pixel, scale, df,ddymin,ddy max, sec, tstr.cnt) 

repeat 

else 

endif 

rent = ent 
avestr = tstr/rent 

open(unit=14,file='maxmin.out , ,status='new') 
if(resp.EQ.1 )then 

write(1 4,51 0)ddxmin,ddxmax, avestr 
510 format(1 x,e1 6.8,/1 x,e1 6.8,/1x,e1 6.8) 
close(1 4) 

else 

write(1 4, 530)ddymin,ddymax, avestr 
530 format(1x, el 6.8,/lx, el 6. 8,/lx, el 6.8) 
close(14) 

endif 


WRITE(9,*)'D0NE' 

READ(9,*) 

BO 
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SUBROUTINE 

DDX(width, height, pixel, scale, df,ddxmin,ddxmax, sec, 

. tstr,cnt) 

real scale, df 
integer*2 width, height 
character^ pixel(8,3000) 
integer 1st, rsf,tof,bof, coif, coll 
real ddxmin,ddxmax,str(1 000) 
real fringe(1 000),ddymin,ddymax,tstr 
character filename*64 
integer i,k,j,n,l,m,sec,cnt 

* calculation of x-fringe locations 

write(9,*)'Calculating x-fringe locations and ddx' 

width = width*8 

do(i=4,8,4) 

do(k=1 ,1000) 
fringe(k) = 0.0 
repeat 


Isf = 0 
rsf = 0 
n = 1 


do(j=1 , width) 

if(pixel(i,j).EQ.'B'.AND.Isf.EQ.O)then 
Isf = j 

elseif(pixel(i,j).EQ.'W'.AND.Isf.NE.O)then 
rsf = j-1 

fringe(n) = (rsf+lsf)/2.0 
n = n+1 
Isf = 0 
rsf = 0 
else 
endif 
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repeat 

* 

* calculation of d/dx 

★ 

n = 2 

100 if(fringe(n+1).NE.0)then 

str(n) - (2.0*df*scale)/(fringe(n+1)-fringe(n-1)) 
if(str(n).LT.ddxmin)ddxmin = str(n) 
if(str(n).GT.ddxmax)ddxmax = str(n) 
n = n+1 
goto 100 
else 
endif 

m = (sec*8-8)+i 


do(k=2,n-1 ) 


write(1 0,500)fringe(k).m,str(k) 
tstr = tstr + str(k) 
cnt = cnt + 1 

repeat 

500 format(1 x,f6.1 ,2x,i5,2x,e1 6.8) 
repeat 


width = width/8 

RETURN 

END 


SUBROUTINE 

DD Y(width, height, pixel, scale, df,ddymin,ddymax, sec, 
. tstr, cnt) 
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integer*2 width.height 
character^ pixel(8,3000) 
integer 1st, rsf,tof,bof, coif, coll 
real df, scale 

real ddxmin,ddxmax,str(1000) 
real fringe(1 000),ddymin,ddymax,tstr 
character filename*64 
integer i,k,j,n,l,m,sec,cnt 

* calculation of y-fringe locations 

write(9,*)'Calculating y-fringe locations and ddy' 

width = width*8 

do(j=4,8,4) 

do (k= 1,1 000) 
fringe(k) = 0.0 
repeat 


tof = 0 
bof = 0 
n » 1 


do(i=1 .height) 

if(pixel{j,i).EQ.'B'.AND.tof.EQ.O)then 
tof = i 

elseif(pixel(j,i).EQ.'W’.AND.tof.NE.O)then 
bof = i-1 

fringe(n) = (bof+tof)/2.0 
n = n+1 
tof = 0 
bof = 0 
else 
endif 
repeat 

* 

* calculation of d/dy 
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n = 2 

110 if(fringe(n+1).NE.0)then 

str(n) = (2.0*drscale)/(fringe(n+1)-fringe(n-1)) 
if(str(n).LT.ddymin)ddymin = str(n) 
if(str(n).GT.ddymax)ddymax - str(n) 
n = n+1 
goto 110 
else 
endif 

m = (sec*8-8)+j 


do(k=2,n-1 ) 


write(1 2,520)m,fringe(k),str(k) 
tstr = tstr + str(k) 
cnt = cnt + 1 

repeat 

520 format(1 x,i5,1 X,f6.1 ,1 X, el 6.8) 
repeat 

width = width/8 
RETURN 

END 


Program STRREG 


implicit none 
integer toolbx 


include quickdraw.inc 
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include file.inc 
include memory. inc 
include misc.inc 
include utilities. inc 
include params.inc 
include grafport.inc 

integer pageheight 
parameter (pageheight=720) 

* An array to contain the bitmap data. 

integer*1 srcbits(72, pageheight) 

* An single line of white pixels (zeros) to fill out the rest of the 

* MacPaint bitmap. 

integer*! srcwht(72) 

integer*1 dstbuf(512) 
integer*1 grafport(grafsize) 
integer grafptr 
integer dstbytes 

* Pointers to the source bitmap and destination buffer. 

integer srcptr, dstptr 
integer i, j, scaniine 

* Bitmap definitions. 

integer bitptr 
integer*1 bitmap(14) 
integer*2 rowbytes 
integer*2 bounds(4) 
equivalence (bitptr, bitmap(l)) 
equivalence (rowbytes, bitmap(5)) 
equivalence (bounds, bitmap(7)) 


character filename*64 
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integer err 

Rectangle for graphics operations on paint file bitmap. 
integer*2 rect(4) 

integer* 1 fraction, reduction, resp.horzshft.vertshft 

* Handle to the VisRgn. 

integer vishandle 

★ 

* Inputs to set region and picture scaling 

★ 

write(9,*)’lnput the region factor, (> 9)’ 

read(9,*)fraction 

write(9,*)" 

write(9,*)'!nput the reduction factor, (< 4)' 

read(9,*)reduction 

write(9,*)" 

write(9,*)’lnput the horizontal shift factor' 

read(9,*)horzshft 

write(9,*)" 

write(9,*)'lnput the vertical shift factor' 

read(9,*)vertshft 

write(9,*)'' 

write(9,*)'Choose the desired mapping' 
write(9,*)'1. Axial Strains’ 
write(9,*)'2. Transverse strains’ 
read(9,*)resp 

write(9,*)" 

write(9,*) 'Now creating the Strain Region MacPaint file.' 

• Initialize the bitmap buffer, grafport. 

do (i = 1, pageheight) 
do (j = 1, 72) 
srcbits(j,i) = 0 
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repeat 

repeat 

* Initialize the line of white pixels. 

do (i = 1 , 72) 
srcwht(i) * 0 
repeat 

* Initialize the grafport record. 

do (i = 1, grafsize) 
grafport(i) = 0 
repeat 

grafptr = toolbx(PTR, grafport) 
srcptr = toolbx(PTR, srcbits) 
cal! too!bx(OPENPORT, grafptr) 

* install a pointer to the bitmap buffer. 

bitptr = srcptr 

* 72 bytes per bitmap row. 

rowbytes = 72 

* Bitmap bounds rectangle. 

bounds(l) = 0 
bounds(2) = 0 
bounds(3) = pageheight 
bounds(4) = 576 

* Make the newly created grafport current. 

call toolbx(SETPORT, grafptr) 


* 


Attach the newly created bitmap to the current grafport. 
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call toolbx(SETPORTBITS, bitmap) 

* Set the PortRect of the new grafport. 

call toolbx(PORTSIZE, 576, pageheight) 

* Set the VisRgn. 

vishandle = long(grafptr + visrgn) 

call toolbx(SETRECTRGN, vishandle, 0, 0, 576, pageheight) 

* Perform some graphics operations on the grafport. 

call STRR EG(f raction, reduction, resp,horzshft,vertshft) 

* The grafport operations are complete. The following code creates 
a 

* MacPaint file from this grafport, and requires toolbx.sub 

* version 2.1 or later. 


* Set the file type and owner. 

do (i = 1, 80) 
params(i) = 0 
repeat 

filename = char(6) // 'strain' 
ionameotr = toolbx(PTR, filename) 

err = toolbx(PBFILECREATE, toolbx(PTR, params)) 
err = toolbx(PBGETFILEINFO, toolbx(PTR, params)) 
fdtype = 'PNTG' 
fdcreator = 'MPNT' 

err = toolbx(PBSETFILEINFO, toolbx(PTR, params)) 

* Open the file. Bitmap data is binary, so use unformatted access. 

open(20, file='strain', status='old', 

+ form='unformatted') 


103 


* Clear out the working buffer for compressing the bitmap data. 

do (i = 1, 512) 
dstbuf(i) = 0 
repeat 

* Write the MacPaint header. 

write(20)(dstbuf(i), i=1 ,51 2) 

* Get a pointer to the bitmap buffer. 

srcptr - toolbx(PTR, srcbits) 

* Compress the data a scan line at a time, writting each compressed 
scan 

* line to the file. 

* 

do (scanline « 1, pageheight) 
dstptr = toolbx(PTR, dstbuf) 
call toolbx(PACKBITS, srcptr, dstptr, 72) 

— UJljJii i.*JUiuX(PTR, dStbuf) 
write(20)(dstbuf(i), i = 1, dstbytes) 
repeat 

* Fill out the rest of the file with compressed lines of white pixels. 

srcptr = toolbx(PTR, srcwht) 
dstptr = toolbx(PTR, dstbuf) 
call toolbx(PACKBITS, srcptr, dstptr, 72) 
dstbytes = dstptr - toolbx(PTR, dstbuf) 
do (scanline = 1, 720-pageheight) 
write(20)(dstbuf(i), i - 1, dstbytes) 
repeat 

close(20) 


end 
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SUBROUTINE STRREG(fraction,reduction,resp,horzshft,vertshft) 

IMPLICIT NONE 

INTEGER xl ,y1 ,x2,y2,N,colm 
INTEGER*4 TOOLBX, size, style, bold, normal, geneva 
INTEGER‘1 black(8),dk1 (8),dk2(8),gray(8),ltgray(8),resp 
INTEGERS Itl (8), It2(8),lt3(8),white(8), fraction, reduction 
INTEGER*1 horzshft.vertshft 

INTEGER*2 rblack(4),rdk1 (4),rdk2(4),rgray(4),rltgray(4) 
INTEGER*2 rltl (4),rlt2(4),rlt3(4),rwhite(4) 

REAL svl ,sv2,sv3,sv4,sv5,sv6,sv7,sv8,sv9,sv1 0 
REAL max,min,temppointy,temppointx,strain1 ,strain2 
character*256 string 


include quickdraw.inc 

black(l) = b'l 1 1 1 1 1 1 1 * 
black(2) = b’l 1 1 11111’ 
black(3) = b*1 1111111* 
black(4) = b'l 1 1 1 1 1 1 1 ’ 
black(5) = b'l 1 1 1 1 1 1 1 ' 
black(6) = b'l 1 1 1 1 1 1 1 ' 
black(7) = b’l 1111111' 
black(8) = b'l 1111111* 

dk2(1) = b'llllllir 
dk 2 ( 2 ) = b'llllllir 

dk2(3) = b'001 10011' 
dk2(4) = b'001 10011' 
dk2(5) = b'l 1 1 1 1 1 1 1' 
dk2(6) = b'l 1 1 1 1 1 1 1' 
dk2(7) = b'001 10011' 
dk2(8) = b'001 10011' 

dkl (1) = b'l 1001 100' 
dkl (2) = b'001 10011' 
dkl (3) = b'l 1001 100’ 
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dkl (4) = 

b'001 10011' 

dkl (5) = 

b'1 1001 100' 

dkl (6) = 

b'001 10011' 

dkl (7) = 

b’1 1001 100' 

dkl (8) = 

b’001 10011' 

gray(1) 

= b'1 1001 100' 

gray (2) 

= b'1 1001 100’ 

gray (3) 

= b'001 10011' 

gray (4) 

= b'001 10011' 

gray(5) 

= b’1 1001 100' 

gray (6) 

= b’1 1001 100' 

gray (7) 

= b'001 10011' 

gray (8) 

= b'001 10011' 

ltgray(1) 

= b’1 010101 O’ 

ltgray(2) 

= b'01010101' 

ltgray(3) 

= b*1 01 01 01 O’ 

ltgray(4) 

= b'01010101' 

ltgray(5) 

= b'1 01 01 01 O' 

ltgray(6) 

= b'01010101’ 

ltgray(7) 

= b'10101010' 

ltgray(8) 

= b'01010101' 

Itl(l) = 

b'1 0001 000' 

Itl (2) = 

b'001 0001 O' 

Itl (3) = 

b'1 0001 000’ 

Itl (4) = 

b'001 00010’ 

Itl (5) = 

b’1 0001 000' 

Itl (6) = 

b'001 0001 O' 

Itl (7) = 

b'1 0001 000’ 

Itl (8) = 

b'001 0001 O' 

It2(1) = 

b'1 0001 000’ 

It2(2) = 

b'OOOOOOOO’ 

It2(3) = 

b'001 00010’ 

It2(4) = 

b’OOOOOOOO’ 

It2(5) = 

b’1 0001 000' 

It2(6) = 

b'OOOOOOOO' 

It2(7) = 

b'001 00010’ 
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It2(8) = b'00000000' 


It3(1) = b'1 0000000' 
It3(2) = b'00000000' 
It3(3) = b'00001000' 
It3(4) = b'00000000' 
It3 (5) = b'1 0000000’ 
It3 (6) = b'00000000’ 
It3(7) = b'00001000' 
It3(8) = b'00000000' 


white(1 ) 
white(2) 
white(3) 
white(4) 
white(5) 
white(6) 
white(7) 
white(8) 


b'1 0000000' 
b’00000000' 
b'00001000' 
b'00000000' 
b'1 0000000' 
b’00000000' 
b'00001000' 
b’00000000’ 


open(unit=1 2,file='maxmin.out',status='old') 

read(1 2,*)min 
read(1 2,*)max 

svl = min 

sv2 = sv1+(max-min)/fraction 
sv3 = sv2+(max-min)/fraction 
sv4 = sv3+(max-min)/fraction 
sv5 = sv4+(max-min)/fraction 
sv6 = sv5+(max-min)/fraction 
sv7 = sv6+(max-min)/fraction 
sv8 = sv7+(max-min)/fraction 
sv9 = sv8+(max-min)/fraction 
svIO = max 


data rlt3 7645,10,655,20/ 
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data rlt2 /660.1 0,670,20/ 
data rltl /675.1 0,685,20/ 
data rltgray /690,1 0,700,20/ 
data rgray /645, 250, 655, 260/ 
data rdkl /660, 250, 670, 260/ 
data rdk2 /675, 250, 685, 260/ 
data rblack /690, 250, 700, 260/ 

data geneva /3/ 
call toolbx(TEXTFONT, geneva) 
data bold /I/ 
style = bold 

call toolbx(TEXTFACE, style) 
if(resp.EQ.1 )then 

string = char(23)//' AXIAL STRAIN BY REGIONS' 
elseif(resp.EQ.2)then 

string = char(28)//'TRANSVERSE STRAIN BY REGIONS' 
else 
endif 

call toolbx(MOVETO,200,10) 

call toolbx(DRAWSTRING, string) 

data normal /0/ 

style = normal 

call toolbx(TEXTFACE, style) 

data size /9/ 

call toolbx(TEXTSIZE.size) 


call toolbx(FILLRECT,rlt3,lt3) 
call toolbx(FRAM ERECT, rlt3) 
call toolbx(FILLRECT,rlt2,lt2) 
call toolbx(FRAM ERECT, rlt2) 
call toolbx(FILLRECT,rlt1 ,lt1) 
call toolbx(FRAM ERECT, rltl) 
call toolbx(FI LLRECT, rltgray, Itg ray) 
call toolbx(FRAMERECT, rltgray) 
call toolbx(FILLRECT, rgray, gray) 
call toolbx(FRAMERECT, rgray) 


108 


call toolbx(FILLRECT,rdk1,dk1) 
call toolbx(FRAMERECT,rdk1) 
call toolbx(FILLRECT,rdk2,dk2) 
call toolbx(FRAM ERECT, rdk2) 
call toolbx(FILLRECT,rblack, black) 
call toolbx(FRAMERECT,rblack) 

call toolbx(MOVETO,30,655) 

WRITE(*,550)sv2,sv3 
call toolbx(MOVETO, 30,670) 

WRITE(*,550)sv3,sv4 
call toolbx(MOVETO,30,685) 

WRITE(*,550)sv4,sv5 
call toolbx(MOVETO,30,700) 

WRITE(*,550)sv5,sv6 
call toolbx(MOVETO,270,655) 

WRITE(*,550)sv6,sv7 

call toolbx(MOVETO,270,670) 

WRITE(*,550)sv7,sv8 

call toolbx(MOVETO,270,685) 

WRITE(*,550)sv8,sv9 

call toolbx(MOVETO,270,700) 

WRITE(*,550)sv9,sv1 0 
550 format(1x, 'strain = ',e10.4,' to ’,e10.4) 

if(resp.EQ.1)then 

open(unit=1 1 ,file='ddy. output', status='old') 
elseif(resp.EQ.2)then 

open(unit=1 1 ,file='ddx. output’, status='old') 
else 
endif 

if(resp.EQ.1)then 

read(1 1 ,*,end=600,ERR=600)x1 ,temppointy,strain1 
yl = ifix(temppointy) 
elseif(resp.EQ.2)then 

read(1 1 ,*,end=600,ERR=600)temppointx,y1 ,strain1 
xl = ifix(temppointx) 
else 
endif 
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xl = ifix(x1 /reduction) - horzshft 

yl = ifix(y1 /reduction) + vertshft 

if(resp.EQ.1)then 

read(1 1 > *,end=600,ERR=600)x2,temppointy,strain2 
y2 = ifix(temppointy) 
elseif(resp.EQ.2)then 

read(1 1 ,*,end=600,ERR=600)temppointx,y2,strain2 
x2 - ifix(temppointx) 
else 
endif 

x2 = ifix(x2/reduction) - horzshft 

y2 = ifix(y2/reduction) + vertshft 


400 if (strain 1 .GE.svl .AND.strainl .LT.sv2)then 
CALL TOOLBX(PENPAT,white) 
elseif(strain1 .GE.sv2. AND.strainl .LT.sv3)then 
CALL TOOLBX(PENPAT,lt3) 
elseif(strain1 .GE.sv3. AND. strainl ,LT.sv4)then 
CALL TOOLBX(PENPAT,lt2) 
elseif(strain1 .GE.sv4. AND.strainl .LT.sv5)then 
CALL TOOLBX(PENPAT,lt1) 
elseif(strain1 .GE.sv5. AND.strainl .LT.sv6)then 
CALL TOOLBX(PENPATJtgray) 
elseif(strain1 .GE.sv6. AND.strainl .LT.sv7)then 
CALL TOOLBX(PENPAT,gray) 
elseif (strainl .GE.sv7. AND.strainl .LT.sv8)then 
CALL TOOLBX(PENPAT,dk1 ) 
elseif(strain1 .GE.sv8. AND. strainl .LT.sv9)then 
CALL TOOLBX(PENPAT,dk2) 
else 

CALL TOOLBX(PENPAT, black) 
endif 

CALL TOOLBX(MOVETO,x1 ,y1) 

CALL TOOLBX(LINETO,x2,y2) 
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xl = x2 

yi = y2 

strain 1 = strain2 
if(resp.EQ.1 )then 

read(1 1 ,*,end=600,ERR=600)x2,temppointy,strain2 
y2 = ifix(temppointy) 
elseif(resp.EQ.2)then 

read(1 1 ,*,end=600,ERR=600)temppointx,y2,strain2 
x2 = ifix(temppointx) 
else 
endif 

x2 = ifix(x2/reduction) - horzshft 
y2 = ifix(y2/reduction) + vertshft 

if(resp.EQ.1)then 

if (x 1 .NE.x2)then 
xl = x2 

yi =y2 

strainl = strain2 
if(resp.EQ.1 )then 

read(1 1 ,*,end=600,ERR=600)x2,temppointy,strain2 
y2 = ifix(temppointy) 
elseif(resp.EQ.2)then 

read(1 1 ,*,end=600,ERR=600)temppointx,y2,strain2 
x2 = ifix(temppointx) 
else 
endif 

x2 = ifix(x2/reduction) - horzshft 
y2 = ifix(y2/reduction) + vertshft 

else 

endif 


elseif(resp.EQ.2)then 


iT^yl ,NE.y2)then 
xl = x2 

yi = y2 

strainl = strain2 
if(resp.EQ.1)then 

read(1 1 ,*,end=600,ERR=600)x2,temppointy,strain2 
y2 = ifix(temppointy) 
elseif(resp.EQ.2)then 

read(1 1 ,*,end=600 ) ERR=600)temppointx,y2,strain2 
x2 = ifix(temppointx) 
else 
endif 

x2 = ifix(x2/reduction) - horzshft 
y2 = ifix(y2/reduction) + vertshft 

else 

endif 

else 

endif 


goto 400 


600 RETURN 
BSD 


